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We have studied the symmetric and non-symmetric pairwise ground state and thermal entangle- 
ment in three-qubit anisotropic Heisenberg (XXZ) and Ising in a magnetic field models and in the 
presence of Dzyaloshinskii-Moriya (DM) interaction. We have found that increasing of the DM inter- 
action and magnetic field can enhance and reduce the entanglement of system. We have shown that 
the non-symmetric pairwise has higher value concurrence and critical temperature (above which the 
entanglement vanishes) than the symmetric pairwise. For negative anisotropy the non-symmetric 
entanglement is a monotonic function of DM interaction while for positive anisotropy it has a max- 
imum versus DM parameter and vanishes for larger values of DM interaction. The conditions for 
the existence of thermal entanglement are discussed in details. The most remarkable result happens 
at zero temperature where 3-qubit ground state entanglement of the system (in spite of 2-qubit 
counterpart) shows the fingerprint of quantum phase transition for an infinite size system. 

PACS numbers: 03.65.Ud, 03.67.Lx, 75.10.Jm, 05.50.-|-q 
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INTRODUCTION 

Entanglement is a property of quantum state which 
has been studied intensively in recent years as a specific 
nonlocal quantum mechanical correlationP, 0, Q and it 
becomes recently a key feature of quantum information 
theory [1]. When the entanglement is generated, obser- 
vation or manipulation of it in practice constitutes a 
major obstacle, because of the fragility of quantum en- 
tanglement to the decoherence induced by environment. 
Therefore, how to generate, maintain and control the en- 
tanglement in the presence of dissipative coupling of the 
system to the environment is of utmost importance in 
the implementation of quantum information processing. 

However, similar to the superposition of two coher- 
ent electromagnetic waves which enable us to learn some 
global information from a localized spatial area, people 
expect that the entanglement, which roots in the same su- 
perposition principle, can enable us to learn some global 
properties from a small part of the system. This obser- 
vation may be one of the main motivations in the recent 
studies d, @, 01 on the role of entanglement between a 
small part, e.g. a block consisting of one or more sites, 
and the rest of the system in the quantum phase transi- 
tion. These results suggested that the local entanglement 
may be used as a good marker of quantum phase transi- 
tion. 

Then over the past few years there has been an on- 
going effort to characterize the entanglement properties 
of condensed matter systems and apply them in quan- 
tum information. The quantum entanglement in solid 
state systems such as spin chains is an important emerg- 
ing fieldji, i, M, ini, 
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Spin chain 

are natural candidates for realization of entanglement, 
and spin effects have been investigated in many other 



systems, such as superconductors [17[, quantum dots|18j 
and trapped ions[l9|. A most known models in the spin 
chains is Heisenberg model and Ising model as a special 
case of Heisenberg model. The Heisenberg model can 
describe interaction of qubits not only in solid physical 
systems but also in many other systems such as quan- 
tum dots fiol nuclear spin 21 1, cavity QED [il. [23l| . opti- 
cal lattice|2J], quantum computation [25| and controUed- 
Not gatejij 

In recent years the two-qubit thermal entanglement 
which includes spin-spin interactions 
27, [2I], and spin-orbit coup ling pol. 
34 1 fDzvaloshinskii(35j-Moriva|36l| interaction) has been 
studied. Entanglement in two-qubit state has been well 
studied in the literature alon g va rious kind of three- 
qubit entanglement states [13, lH, [H, \M, El- The 
three-qubit entanglement states have been shown to 
possess advantage over the two-qubit states in quan- 
tum teleportation[4^ . dense coding[4^ and quantum 
cloning[44|. More specificly, both quantitative and qual- 
itative behavior of the entanglement in two-qubit and 
three-qubit systems are different. The concurrence is al- 
ways increasing in a two-qubit model versus DM param- 
eter while it is both increasing and decreasing versus DM 
for a three-qubit model. The position where the concur- 
rence is maximum or zero for a three-qubit model corre- 
sponds to the quantum critical point of the infinite size 
system while it is not the case for a two-qubit model. 
Moreover, the pairwise entanglement can be defined as 
symmetric and non-symmetric ones with differente prop- 
erties for three-qubit as will be discussed in this article. 

In addition to the above facts, recently some novel 
magnetic systems with antiferromagnetic (AF) prop- 
erties, such as Cu{C6D5COO)2^D20^, YbAAsa^, 
BaCu2Si207^, a - Fe203, LaMnOa^ and 
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-K'2V308[49|, were discovered in the category of 
quasi-one dimensional materials which are known 
to belong to an antisymmetric interaction of the form 
D .{Si X Sj) which is known as the Dzyaloshinskii- 
Moriya (DM)interaction[3^, [3^1 ■ Thus, investigation of 
the quantum effects in such spin models requires more 
research in this direction. 

In this paper, we have investigated the influence of 
the the anisotropy coupling, magnetic field and the z- 
component DM interaction on the non-symmetric (pi2 = 
Tr^piT)) and symmetric pairwise (pi3 = Tr2p{T)) 
ground state and thermal entanglement of the three- 
qubit XX Z and Ising models in a magnetic field and 
in the presence of DM interaction. We will show that the 
DM interaction, anisotropy and magnetic field parame- 
ters are efficient control parameters of entanglement. In- 
creasing the DM coupling and anisotropy have a different 
effects on the entanglement and can enhance or reduce 
the entanglement, whereas these parameter just increase 
the entanglement in two-qubit [H, 31, 32, 33, [s^l 
counterpart. However we show that at T = the 3-qubit 
quantum phase transition points of this system corre- 
spond to those ones at thermodynamic limit 15, 50j . 



THE MODEL INTERACTION 

The Hamiltonian of N-qubit XXZ model with z- 
component of DM Interaction is 

~ N 

H{.l A) = - ^ [afaf+i + + ^crlat^, (1) 

i 

where J is the exchange coupling, D is the strength of 
z component DM interaction and A defines the easy- 
axis anisotropy which can be positive or negative. The 
positive or negative J corresponds to the antiferromag- 
netic (AF) or ferromagnetic (F) cases, respectively, erf 
refers to the a-component of Pauli matrix at site i. A 
TT-rotation around z axis on odd (or even) sites maps the 
F case (J < 0) to the AF case with the opposite sign of 
anisotropy, 

J N 

H{J, A) = - ^ [af af+i + af+, - Aafa^^, (2) 

i 

+Dia^al,-afa:^,)\, J ^ \J\ > 0. 

So we can restrict ourselves to AF case (J > 0) with 
D > and arbitrary anisotropy (A < and A > 0) 
without loss of generality. 

We can also restore the Ising model with DM interac- 
tion (IDM) in the limit J — > 0, A — * oo and D oo 



where J A - 
is given by 



J and ^ 



D. The resulting Hamiltonian 



i N 



i) • (3) 



It should be noticed that both AF (corresponds to the 
positive anisotropy in XXZ with DM) and F cases (cor- 
responds to the negative anisotropy in XXZ with DM) 
have to be considered separately, since the mentioned 
symmetry can not be justified for Eq.Q. We will con- 
sider the above Hamiltonians (Eqs.([2l [3l)) on a string of 
3-qubits and obtain the entanglement of them in terms 
of the model parameters. 

The entanglement of two-qubit can be measured by the 
concurrence which is defined aslSlI 



Cij = max |2 max(AA;) 
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fc=i 



.,0} 



(4) 



where Afc (k=l,2,3,4) are the square roots of the eigen- 
values of the following operator 



(5) 



Pij is the density matrix of pair i and j spins and as- 
terisk denotes the complex conjugate. For a system in 
equilibrium at temperature T the state of a system is 
determined by the density matrix 



piT) 



-/3H 



(6) 



where H is the system Hamiltonian, Z — Tr{e~^^) is the 
partition function and (3 = where T is temperature 
and Kb is Boltzmann constant. For simplicity we take 
Kb = 1. 

To get the concurrence of two qubits in a string of 
3-qubits we define two types of reduced density matrix, 
the symmetric and non-symmetric ones. Let label the 
3-qubits as 1,2,3 sequentially. The symmetric reduced 
density matrix (pia) is defined as pi3 = tr2{p), where 
p is the density matrix of 3-qubits. The non-symmetric 
reduced density matrix is pi2 = tr^^p). 



THREE QUBITS XXZ WITH DM INTERACTION 

A straightforward calculation gives the eigenstates and 
eigenvalues of the 3-qubit XXZ-I-DM Hamiltonian [isj . 
The positive square roots of the eigenvalues of the matrix 
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FIG. 1: (color online) The concurrence of non-symmetric pair- 
wise is plotted versus A for different values of T, {D — 0.0). 



FIG. 3: (color online) The variation of non-symmetric pair- 
wise concurrence with DM coupling (D) for different values 
of T, A = -1.5. 




FIG. 2: (color online) The concurrence of non-symmetric 
pairwise against anisotropy (A) for different values of T 
,{D = 2.0). 



The partition function is then 



Z = 4e*^^cosh(^/3(7)+4e-T^^cosh(^/3A) (8) 

The square root of matrix R is invariant under the sub- 
stitution D —D, so the generality of AF case with 
D > and arbitrary anisotropy, remains pristine. 

At T = 0, for A > — Vl + the concurrence of the 
doubly-degenerate ground state is 



,J 



2{q + A- Vl + D^)Vl + 



9(9 + A) 

while it is a disentangled ferromagnetic states for A < 
-Vl + 8D2. At the level crossing line (A = + SD^) 
the ground state is 4-fold degenerate and disentangled. 

For nonzero temperature (T 7^ 0) the concurrence is 



R for non-symmetric pairwise is obtained 



Ai 



1 rl 



z 12 



e 2 



il3A 



T/3A 



i cosh(^/3g) - ^ sinh(^/3g) |, 



A - SVl + . , ,J 
— sm. 



^cosh(^/3g) 

-{- 

Z\2 

- cosh(-/3g) + sm. 



2q 



(7) 



CiY(r)-^max{( 



1 



)smh(-/^g) 

q 4 



- cosh(^/3g) 



2e- 



The model is entangled if T > rc(A, D), and rc(A, D) 
is determined by the following nonlinear equation 



^ l_+_87^^ sinh(^/3g) - cosh(^/3g) = e 



-4/3. 



(9) 



To determine whether the entanglement exists or not, 
we have to consider two different cases. 

• A < 

The ground state entanglement of non-symmetric 
pairwise have been shown in Fig.((T]) and Fig.(l2|) 
as a function of A for D = {] and D = 2 which 
are denoted by T = plots. Both figures mani- 
fest that the concurrence of ground state decreases 
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FIG. 4: (color online) The variations of non-symmetric pair- 
wise concurrence as a function of DM coupling (D) for differ- 
ent values of T, A = —0.5. 



FIG. 5: (color online) The concurrence of non-symmetric pair- 
wise is plotted versus (Z>) for different values of temperature, 
A = 0.5. 



by reduction of the anisotropy and suddenly be- 
come zero below the critical line Ac = — \/l + 1?^. 
In the other words, for each value of £>, there is a 
threshold Ac = — vl + Z?^ above which the ground 
state will be entangled. Accordingly, for each value 
of anisotropy (A < —1) there is a value of DM 
coupling Dc = %/ A^ — 1 under which the ground 
state looses its entanglement. The critical line 
Ac = — VT+D^ is at the position where the level 
crossing occurs and corresponds to the critical line 
of this model at thermodynamic limit (i.e. infinite 
number of qubits). At this critical line the global 
U{1) X Z2 symmetry of the Hamiltonian changes to 
the local SU{2) symmetry. 

At finite temperature, thermal entanglement be- 
haves similar to the ground state counterpart ex- 
cept that it becomes disentangled gradually by in- 
creasing temperature. For T 7^ 0, the critical line 
below which the thermal entanglement vanishes is 
a function of D and temperature, i.e Ac(£',T) 
and is given by Eq.(l9]). Generally, the decrease 
in anisotropy parameter and increase in temper- 
ature have a reduction infiuence on the concur- 
rence. The effect of temperature is plotted in 
Fig.dD) and Fig. ([2]) for some fixed values of tem- 
perature. To scan the influence of DM coupling the 
concurrence has been plotted in Fig. ([3]) and Fig.(j4|) 
at flxed values of A = —1.5 and A = —0.5 versus 
D. For A < — 1, the concurrence jumps suddenly 
to non-zero value as D crosses the critical value 
Dc = VA^ — 1 as shown in Fig. ([3]). Contrary to 
the previous case, for —1 < A < the ground state 
is always entangled while the thermal entanglement 
becomes zero for \D\ < Dc(A, T) which is justifled 
in Fig.©. 




FIG. 6: (color online) The variations of non-symmetric pair- 
wise concurrence with DM coupling (D) for different values 
of T, A = 1.5. 



At A = the concurrence of ground state is inde- 
pendent of DM coupling and take a constant value. 



Ci2(T = 0,A = 0) = 



2V2-I 



• A > 

In the positive anisotropic region, increasing the 
anisotropy enhances the ground state and thermal 
entanglement up to a maximum and then decreases 
gradually which can be seen in Fig. ([5]) and Fig.®. 
At T = the maximum value of concurrence is 
due to the maximum fluctuations which exist at the 
critical line Ac = vl + 8D^. The critical line be- 
longs to the spin-fluid to Neel phase transition for 
infinite number of qubits [l^. Contrary to what 
happened in the negative anisotropy region where 
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the ground state entanglement vanishes immedi- 
ately below the critical line Ac = —y/l + D^, in 
the positive anisotropy region the ground state en- 
tanglement has non-zero value in the both side of 
critical line Ac = Vl + D'^- However, in the both 
cases the global U{1) x Z2 symmetry of the Hamil- 
tonian breaks to the local SU{2) symmetry on the 
critical line (Ac = ^Vl + D^). The influence of 
temperature is to decrease the concurrence and in- 
crease the critical anisotropy parameter {Ac{D, T)) 
beyond which the system being entangled (as given 
by Eq.®). To survey the effect of DM interac- 
tion the ground state and thermal concurrence have 
been plotted in Fig. ([5]) and Fig.® versus D for 
A = 0.5, A = 1.5 and different values of temper- 
ature. As far as < A < 1 the influence of DM 
interaction is to reduce the concurrence of ground 
state (Fig. ([5])) while for A > 1 the ground state 
concurrence gets a maximum at the position of the 
critical line Ac = yl + D^. In the other words, for 
A > 1 increasing of the DM interaction raises the 
ground state concurrence until its maximum value 
and then decreases slowly (Fig.®). This is a fea- 
ture which is only observable in 3-qubits system 
and is related to the fact that the 3-qubits system 
show the critical behavior of infinite qubits limit 
correctly. For low temperatures, thermal entangle- 
ment behaves similar to the ground state concur- 
rence, but at high temperatures the thermal en- 
tanglement is increased on the onset of the DM 
coupling. 

• Symmetric pairwise 
We have discussed extensively the non-symmetric 
pairwise entanglement (C12) which has been calcu- 
lated from pi2. Now, we summerize the main fea- 
tures of the symmetric pairwise entanglement (C13) 
in comparison with the non-symmetric one without 
presenting the details. 

In the negative anisotropy region (A < 0), there is 
a threshold anisotropy and D parameters beyond 
that the symmetric and non-symmetric pairwise 
are entangled. Decreasing the anisotropy increases 
the ground state entanglement of symmetric pair- 
wise but decreases the ground state entanglement 
of non-symmetric pairwise. However, reduction of 
the anisotropy increases the thermal entanglement 
of symmetric pairwise up to a maximum value and 
then decreases it gradually to be vanished while it 
decreases the non-symmetric pairwise thermal con- 
currence to zero. In symmetric pairwise the incre- 
ment of DM coupling decreases the ground state 
entanglement but enhances the thermal entangle- 
ment up to the maximum value and then decreases 
it. For non-symmetric pairwise the influence of DM 
interaction on ground state and thermal concur- 



rence is incremental. For the positive anisotropy 
region (A > 0), the symmetric case is disentangled 
in the presence or absence of temperature. 



THREE QUBITS ISING WITH DM 
INTERACTION 

We have studied both symmetric and non-symmetric 
pairwise entanglement of the Ising model with DM inter- 
action defined by Hamiltonian in Eq.Q. Its qualitative 
behavior is similar to the XXZ model with DM interac- 
tion while only some quantitative changes observed, for 
example the critical point for both AF and F cases is 
_Dc = 1. Hence, we do not present the results of 3-qubits 
Ising model with DM interaction here; however, we will 
show the effect of magnetic field on the entanglement in 
this model. Let rewrite the Hamiltonian of 3-qubit Ising 
model with DM interaction in the following form 



H 



(10) 



where h is proportional to the strength of magnetic field, 
J and D arc the exchange and DM couplings, respec- 
tively. The square root of the eigenvalues of matrix R 
(Eq.®, are invariant under the substitutions D —D 
and h ~h, so we will consider only D > and h > 
without loss of generality. 

As mentioned before we consider the AF and F cases 
of Ising model with DM interaction separately. However, 
we only study the non-symmetric pairwise entanglement 
(C12) because the symmetric one (C13) shows the same 
results qualitatively. 

• AF case (J > 0) 
For zero temperature (T = 0), the ground state is 
an entangled one and its concurrence is given by 



C 



0,/i ^ 0) 



2D 



which is plotted in Fig.® and denoted by T = 
0. The thermal concurrence is nonzero for T < 
Tc{D, h) which is shown for different temperatures 
and h = 2.0 versus D in Fig.®. We have ob- 
served that the the critical temperature decreases 
with increasing of the magnetic field. Moreover, the 
DM interaction creates entanglement in the system. 
For zero temperature, the onset of DM interaction 
leads to nonzero concurrence while for finite tem- 
perature D should be greater than a critical value 
to have nonzero concurrence. The concurrence is 
plotted versus h for different values of temperature 
and D = 2.0 in Fig.®. This figure shows that 
for low temperature the concurrence increases with 



D 



D 



FIG. 7: (color online) The variations of non-symmetric pair- FIG- 9: (color online) The concurrence of non-symmetric pair- 
wise concurrence versus D for different values of T in AF case wise against D for different values of T for F case ( J = -1), 
(J = 1), /i = 2.0. ft = 2.0. 
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FIG. 8: (color online) The concurrence of non-symmetric pair- 
wise as a function of magnetic field (h) for different values of 
T in AF case (J = 1), D = 2.0. 



2D 



For non-zero temperature and D < Dc{T, h) the 
concurrence is zero as shown in Fig. ([9]) for h = 2.0. 
The increment of DM interaction and the magnetic 
field induce entanglement in this system as far as 
T < Tc{D,h). The critical value of DM couphng 
{Dc{T, h)) is increased with the increment of mag- 
netic field. Moreover, we have observed that the 
amount of concurrence is enhanced for higher mag- 
netic fields. More importantly, all figures reveal 
that the system can be entangled in a region which 
is not entangled at T = by the effect of thermal 
fluctuations. 



SUMMERY AND CONCLUSIONS 



the increase of magnetic field initially to reach a 
maximum and then decreases gradually for larger 
magnetic field, whereas for mid-range temperature 
the magnetic field reduces the concurrence slowly. 
In other words, for low temperatures field induces 
entanglement in a 3-qubit system of Ising model 
with DM interaction. 

• F case (J < 0) 
For D < ^[{3 + 2hY - l]/8, the entanglement of 
the ground state is zero while it is nonzero for D > 
^ [(3 -|- 2/i)2 — l] /8. In the entangled region the 



The ground state and thermal entanglement of sym- 
metric and non-symmetric pairwise in three-qubit XX Z 
and Ising models in the presence of DM interaction has 
been investigated, we have studied the influence of DM 
and anisotropy parameters on the concurrence of these 
models. 

The DM interaction and anisotropy are the efficient 
control parameters of entanglement. For the positive 
anisotropy region (AF case of IDM), the symmetric 
pairwise is disentangled in finite and zero temperature. 
For non-symmetric pairwise, increasing of anisotropy en- 
hances the ground state and thermal entanglement of 
the system to its maximum value and then decreases it 
slowly. In this region the DM interaction has a different 
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FIG. 10: (color online) The concurrence of non-symmetric 
pairwise is plotted versus h for different values of T in F case 
(J= -1), D = 2.0. 



effect on the concurrence. For < A < 1, the increase 
of DM coupling decreases the ground state entanglement 
and concurrence at low temperature but at high temper- 
ature DM interaction enhances thermal concurrence. For 
A > 1, increment of DM coupling enhances the ground 
state entanglement and low temperature concurrence to 
a maximum value and then decreases it gradually but it 
has only an increment effect on the concurrence at hight 
temperature. 

In the negative anisotropy region (F case of IDM) 
there is a threshold for anisotropy and D parameters 
above which both symmetric and non-symmetric pairwise 
are entangled. Decreasing the anisotropy increases the 
ground state entanglement of symmetric pairwise but de- 
creases the ground state entanglement of non-symmetric 
pairwise. However, lowering the anisotropy increases 
the thermal entanglement of symmetric pairwise up to a 
maximum value and then decreases gradually to be van- 
ished, but it decreases the non-symmetric pairwise ther- 
mal concurrence to be eliminated. In symmetric pairwise 
the increment of DM coupling decreases the ground state 
entanglement but enhances the thermal entanglement up 
to the maximum value and then decreases it. For non- 
symmetric pairwise the influence of DM interaction on 
the ground state and thermal concurrence is incremen- 
tal. 

The most noteworthy result occurs at T = where 
the 3-qubit ground state entanglement of the systems 
shows the fingerprint of quantum phase transition for 
an infinite size system. For negative anisotropy case of 
both symmetric and non symmetric pairwise, level cross- 
ing of the ground state and first excited state occurs at 
A = - y/{l + D^) {D=l for F case of IDM) under which 
the concurrence is zero and jumps to a nonzero value 



for A > — -I- D^) {D < 1). For positive anisotropy 
case of non-symmetric pairwise, the symmetry break- 
ing of the ground state occurs similar to the negative 
anisotropy case without level crossing at A = ^(1 + D^) 
(D=l for AF case of IDM) and the concurrence has a 
maximum value due to the maximum quantum fluctu- 
ations. At these points the global U{1) x Z2 symme- 
try of the Hamiltonian is changed to local hidden SU (2) 



symmetry [52|, [5 



We have probed the influence of magnetic field in 3- 
qubit IDM model. In non-symmetric pairwise and at low 
temperature the magnetic field induces the entanglement 
to the system which get a maximum for finite field. Ob- 
viously for high magnetic field the model enters a para- 
magnetic phase which is disentangled. For middle range 
of temperature the entanglement is reduced by adding 
the magnetic field, whereas for high temperature, there 
are many states populated and its reduction is very tiny. 
It is related to the thermal mixing of states which is a 
source of entanglement in the system. For the ferromag- 
netic case, the variations of entanglement is faster than 
the AF one. 

For symmetric pairwise in the AF case the magnetic 
field can induce entanglement in a non-entangled system 
{h = 0) and the entanglement enhances with increasing of 
the DM interaction. The entanglement is increased with 
magnetic field to reach a maximum value while further 
increment of magnetic field decreases the entanglement 
which gradually disappears. This happens very quickly 
for F case and the critical temperature where the en- 
tanglement is vanishing is lower than the AF case. The 
optimal mixing of all eigenstates in the system leads the 
maximum value. 

More generally, the entanglement properties of a finite 
system depends on the number of qubits as discussed in 
the text. 

The authors would like to acknowledge M. Kargarian 
and M. F. Miri for useful discussions and comments. This 
work was supported in part by the Center of Excellence in 
Complex Systems and Condensed Matter (www.cscm.ir). 
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